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We study the density of states in graphene at high magnetic field, when the physics is dominated 
by strong correlations between electrons. In particular we use the method of Haldane pseudopo- 
tentials to focus on almost empty or almost filled Landau levels. We find that, besides the usual 
Landau level peaks, additional peaks ("sashes") appear in the spectrum. The energies of these 
peaks are determined by the strength of Haldane's pseudopotentials, but as opposed to the usual 
two-dimensional gas, when there is a one-to-one correspondence between a Haldane pseudopotential 
and a peak in the spectrum, the energy of each peak is determined in general by a combination of 
more than one pseudopotential values. An eventual measure of these peak in the density of states 
spectrum of graphene would allow one to determine the value of the pseudopotentials in graphene, 
and thus test the strength of the interactions in this system. 

PACS numbers: 73.22.Pr, 68.37.Ef, 71.70.Di 



I. INTRODUCTION 

As high quality samples of graphene are now readily 
available^— an exciting prospect is to observe signatures 
of strong electron correlations in this system. Graphene's 
distinctive feature with respect to a more conventional 
two-dimensional (2D) electron gas such as the one re- 
alized in GaAs-based structures is that its charge car- 
riers are relativistic in character, having a linear energy 
dispersions^ From a theoretical perspective it is interest- 
ing to understand the interplay between the Dirac charac- 
ter of the carriers and the effects of interactions between 
the electrons^— While the linearity of the spectrum has 
already been demonstrated by different experiments^ 5 - the 
effects of interactions have been believed to be quite 
small in standard graphene samples. Recently, however, 
the fractional quantum Hall effect has been observed in 
suspended graphene^^ demonstrating that interactions 
can play an important role in graphene in a strong mag- 
netic field with sufficiently high mobility. 

The analysis of fractional quantum Hall (FQH) states 
arising at high perpendicular magnetic field in a two- 
dimensional system has revealed most prominently the 
effects of electronic interactions^ - — In GaAs-based 2D 
electron gases in the FQH regime important informa- 
tion has been obtained from transport and shot noise 
experiments! 14 ' 15 However, it is known that various 
2D correlated states do not show anomalous transport 
signatures ^ Moreover, due to the manner of construc- 
tion of the quantum well structures, having a direct ac- 
cess to their bulk properties is in general not straightfor- 
ward. For example, only a few measurements of the local 
density of states have been performed in 2D gasesr^— , 
in general via indirect methods. 

Compared to the other 2D gases, graphene has the 
advantage that being a surface electronic system, in ad- 
dition to standard transport measurements, it is directly 



accessible by local density of states (DOS) measurements 
such as scanning tunneling microscopy (STM)^S This 
opens the perspective to use such measurements to ob- 
tain information about the electronic interactions. Here 
we study the DOS in graphene exposed to a strong mag- 
netic field focusing on the effects of the electronic inter- 
actions, taken into account via the method of Haldane 
pseudopotentials^ 

Our analysis has been motivated by a recent experi- 
mental analysis of the high-resolution time-domain ca- 
pacitance spectroscopy for a 2D quantum well in the 
quantum Hall regime^ The observation of unexpected 
peaks in the high-energy spectra for LL filling factors 
near v = and v = 1, to which the authors referred as 
"sashes" , has called for new perspectives in the problem 
of 2D correlation physics in strong magnetic fields. The 
filling factor v = n e i/ns is the ratio between the elec- 
tronic density n e i and that ns = eB/h of flux quanta 
threading the 2D surface of the system. Soon after the ex- 
perimental observation, a few different approaches have 
been proposed for the description of these peaks in rela- 
tion to the electronic interactions . 21 ' 22 In particular, in 
Ref. Hi] these peaks have been attributed to the strong 
correlations between electrons, which have been modeled 
using the Haldane's pseudopotentials. 

In the present paper, we generalize this calculation for 
2D electrons in graphene. In this case, the total filling 
factor is defined as v = v n + v, where v n = ±2(2?? + 1) = 
±2, ±6, ... is the filling at which relativistic integer quan- 
tum Hall effect occurs^ and v is the partial filling factor 
of the n-th LL. In this approach, at very low partial filling 
factors v — > 0, besides a number of completely filled and 
inert LLs, only one extra particle is considered present in 
the ground state. Whereas this theoretical limit allows 
for an important simplification in the calculation of the 
tunneling density of states, it is experimentally relevant 
as long as the average distance between particles d ~ 
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Ib/Vv of electrons in the n-th LL is larger than the cy- 
clotron radius Rc = Ib \/^n + Snfi in graphene, in terms 
of the magnetic length lg = ^JficjeB ~ 26 ran/ y/ B[T]. 
Our theoretical analysis is therefore applicable in the 
limit v < l/(27i + <y nj0 ). 

Assuming thus that the state of the system is a one- 
particle state, when a second particle tunnels in from the 
STM tip, the measured STM signal is proportional to 
the overlap between the resulting state, which is not an 
cigenstatc of the high-magnetic field two-particle Ham.il- 
tonian, and the two-particle eigenfunctions. These eigen- 
functions are characterized by a quantum number asso- 
ciated with the two-particle relative angular momentum. 
As a result, discrete peaks arise in the density of states^ 
corresponding to the energy differences between the two- 
body interacting eigenfunctions of the system and the 
one-particle state. These energy differences are given by 
Haldanc's pseudopotcntials, and their measurement via 
the peaks in the spectrum yields information about the 
interaction parameters in the system. 

Similarly to non-relativistic 2D electron systems, such 
as in GaAs heterostructures, we find that the high-field 
DOS, close to the filling v n = ±2(2n + 1), allows for a 
determination of Haldanc's pseudopotcntials and thus of 
the effective Coulomb interaction in graphene. However, 
contrary to non-relativistic 2D electrons in a strong mag- 
netic field, the center-of-mass (CM) is not separable from 
the relative coordinate as a consequence of the Lorentz 
invariance. The extraction of Haldanc's pseudopotentials 
in higher LLs (n ^ 0) is therefore more involved than in 
n = 0, where the two-body problem is equivalent to that 
in non-relativistic 2D electron systems. Our results may 
be tested in high-field STM that has already been ap- 
plied successfully in the past to graphene^, and which 
has also been proposed as a tool for the study of high-field 
electron-solid phases^, as well as of the role of impurity 
scattering £2r2IL 

The paper is organized as follows: in Sec. UH we re- 
view the recent theoretical interpretation of the "sash" 
features observed in high-quality 2DEG samples via Hal- 
dane's pseudopotentials. In Sec. IIIIl we generalize this 
theory to graphene by solving the two-body problem for 
Dirac quasiparticles in a strong magnetic field. With the 
help of the exact two-body eigenstates of the interacting 
system, in Sees. IIVI and fVl we calculate the DOS, and we 
describe how it can be used to extract information about 
the pseudopotentials for the LLs n — and ti/0, respec- 
tively The last section (Sec. IVI|) of the paper presents a 
discussion of the results and the conclusions. 



II. THEORETICAL INTERPRETATION OF 
THE DOS UNDER HIGH MAGNETIC FIELD 
USING HALDANE PSEUDOPOTENTIALS 

Here we review the main aspects of the theory based 
on Haldanc's pseudopotcntials proposed in Ref. [2l| to ex- 
plain the unexpected sashes that have been recently ob- 



served in the DOS of a 2D gas at high magnetic fields 
In standard STM experiments, the measured differential 
conductance is taken to be proportional to the tunneling 
DOS of the system being probed. At zero temperature, 
this quantity is given by 

A(w) = J2 K*a(# + ^)\4^(N))\ 2 S(uj - E afi ) 

Oi 

+ Y,\(V a (N - l)\c p \y (N))\ 2 5(Lo + E afi ), (1) 



where |\I , o(A r )) is the ground state of the iV-particle sys- 
tem and c~P is a fermionic operator that removes (adds) a 
particle from (to) the one-particle state labeled by a set of 
quantum numbers /3. In the high-field case for graphene, 
which we consider here, the set of quantum numbers is 
given by f3 = {n, m, a, K/K'}, where n is the LL index, m 
labels the degenerate one-particle states inside the level, 
a is the spin index and K/K' is the valley index^i Thus, 
the first term describes an electron added to the system, 
whereas the second one corresponds to removing an elec- 
tron. Here, \^> a (N ± 1) ) are the exact eigenstates of the 
interacting N ± 1-particle system labeled by another set 
of quantum numbers a. Their energy difference is given 
by E a _ = [E a (N ± 1) - E (N)}. Because the quan- 
tum number m is associated with the center of the elec- 
tronic cyclotron motion, which is a constant of motion, 
the DOS is independent of m in a translationally invari- 
ant system<2i Since the DOS yields information about 
the spectral properties of the many-body system, STM 
experiments can thus provide important information on 
the correlation physics in the quantum Hall regime. 

To evaluate the above formula, in Ref. Hi], one consid- 
ers the extremely dilute limit (y ~ 0), where the ground 
state simply consists of a single particle in the lowest LL 
|\f'o(A r = 1)). Due to particle-hole symmetry, these re- 
sults also hold for the electron-removal part of the DOS 
when the system is close to an almost filled LL state. The 
tunneling experiment in this regime is then described by 
adding an extra particle instantaneously to the ground 
state (without perturbing it). The resulting state is not 
an eigenstate of the Hamiltonian, the exact two-body 
eigenstates \^ a (N = 2)) of a 2D gas in a strong mag- 
netic field being labeled by their relative angular mo- 
mentum. Moreover, the eigenvalues corresponded to the 
two-particle interacting states are shifted by the Coulomb 
interactions, the corresponded shifts being described by 
the Haldane's pseudopotentials^. The resulting spectral 
peaks in the DOS spectrum occur at energies matching 
the difference between the energy of the interacting two- 
particle states and that of of the one-particle states - this 
difference being precisely given by the Haldane's pseu- 
dopotentials. Thus the DOS spectrum provides detailed 
information about the values of the pseudopotentials in a 
system, which is crucial for understanding the formation 
of various correlated ground states. 
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III. TWO-BODY EIGENSTATES OF DIRAC 
PARTICLES UNDER MAGNETIC FIELD 

As discussed in the previous section, knowledge of the 
exact two-body interacting eigenstates and their corre- 
sponding eigenvalues, as well as of the non-interacting 
state resulting by the addition of an extra particle to 
the one-particle state, is required in order to calculate 
the DOS in the low-filling-factor limit. Given that the 
two-particle eigenstates of the non-interacting Hamilto- 
nian are also eigenstates of the Coulomb interaction po- 
tential, the two-particle interacting eigenstates have the 
same form as the non-interacting ones, but correspond to 
different eigenvalues. They are slightly more complicated 
for graphene for which the relevant charge excitations are 
Dirac-like particles, than for a 2D gas. Besides, one needs 
to take into account the extra spin and valley degrees of 
freedom. To begin, in Sect. IIIA, we focus on the orbital 
part of the interacting wavcfunction for two Dirac par- 
ticles in a magnetic field. In Sect. IIIB we write down 
the wavefunction obtained when adding an extra Dirac 
particle to a one-particle state. The effects of the extra 
degrees of freedom (spin, valley) will be touched upon in 
Sect. IIIC 



A. Two-particle interacting eigenstates and the 
corresponding eigenvalues 

As mentioned above the two-particle interacting eigen- 
states are the same as the non-interacting ones. In order 
to calculate them, one notes that at low-energy (up to 
fractions of an eV) , the electronic properties of graphene 
can be described using a continuum model, with elec- 
trons localized around the two Dirac cones, convention- 
ally called the K and K' valleys £ The Schrodinger equa- 
tion for a Dirac particle in the K valley in a magnetic 
field can be written as 



Hm = va- = 



(2) 



where the Fermi velocity v ~ 10 6 m/s plays the role of 
the speed of light. Furthermore, a = (a x 1 a v ) in terms 
of Pauli matrices, and II = (II^, 11^) is the canonical mo- 
mentum operator, after Peierls substitution, which obeys 
[n^jlly] = —ih 2 /l 2 B . By defining the ladder operators 
a = l B (n x - m y )/V2h and a) = l B (J^ x + iIL)/y/2h 
that obey the harmonic oscillator algebra, [a, a'] = 1, 
the Schrodinger equation becomes 



V2h 



a 
^ 



* = E^. 



(3) 



Diagonalizing the Hamiltonian yields the energy spec- 
trum of rclativistic LLs 



E 



\h—V2n, 
Ib 



(4) 



where A = ±, and n = 0,1,... denotes the LL. The 
associated eigenstates are given by 





n = 0, m] 



for n = 



1 / \n — l,m) 
~ V2\ A|n,m) 



*An,m = ( ''"il^r I , for 71 ^ 0, (5) 



in terms of the (non-relativistic) LL states |n, m) = 

■^-^==— 10). Here, we have implicitly introduced the 

cyclotron-orbit-center operator b = [x + iy H — — (H x + 
iTly)] that has the associated quantum number m with 



tfb I ®r. 



m^n^n, where I is the 2x2 identity ma- 



trix. Since [a, b] = [a\b] = 0, and [H,tfb I] = 0, the 
additional quantum number m labels the macroscopic 
degeneracy Nb = n B A for each Landau level n, in terms 
of the total surface A. 

The solutions of the same problem for the second val- 
ley K' can be obtained by the transformation ^ — > a x ^ , 
since the the Hamiltonian around the K' valley is related 
to Eq. (2) via H — > a x H<r x . This simple covariance prop- 
erty permits to focus the discussion in the rest of this only 
on one K valley (assuming that the interaction is invari- 
ant under the above transformation, for more details see 
Sect. IIIC). 

For the two-body problem of Dirac particles in the ab- 
sence of magnetic field, it was shown^S that due to the 
coupling of sublattice and orbital degrees of freedom, the 
CM degree of freedom cannot in general be separated 
from the relative coordinates degree of freedom. In the 
presence of a magnetic field, we have 



u[o--rii ®i- 

or explicitly 



I <g) a ■ ri 2 ] * = E ty, 



\~~1b~ 



/0 a 2 oi 0\ 



at 
a 2 



(G) 



(7) 



\ a\ 4 / 



where the index 1, 2 labels the two particles. By defining 
the CM and relative coordinates, zr = (z\ + z 2 )/2 and 
z r = z\ — z 2 , respectively, the ladder operators become 
a R = ( a i +fl2)/v / 2, a r = (at — a 2 )/\/2. The Schrodinger 
equation then takes the form 



"r 

a]. 



Cifl — d r 






V o 



a R + a r 



(4 - at 





cir + a r 
o-r — a r 




(8) 



which shows that this is also the case in the presence of 
a magnetic field. 
We denote 



VNW.Mlml 



(9) 
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where we have defined bjj = (bi + 6 2 )/\/2, K = (bi — 
6 2 )/\/2 as the LL ladder operators, with the subscripts R 
and r in \N, M)n\n,m) r indicating the CM and relative 
coordinates subspaces. Using this notation one may write 
down the cigenstatcs in a 4-spinor form \^>M. m (N = 2)), 
see Table U for the exact form of the lowest energy eigen- 
states. We remark that the macroscopic degeneracies in 
both the CM and relative angular momenta, M and m, 
are similar to the single-particle LL case. 

While the eigenstates for the non-interacting and in- 
teracting problems are the same, the values of the eigen- 
values arc different. In order to compute these values, we 
now take into account the interactions 

[w fir + va- ff 2 + V , i l2 I<8>l]* = E^f, (10) 

where we consider the interaction potential to be 
isotropic (n, r 2 \Vi,2 |n, ^2) = V(\n — r 2 |). The eigen- 
values E of the fully interacting Hamiltonian in Eq. (10) 
can be obtained by sandwiching it between the eigen- 
states \^M,m(N = 2) ) described above, with the help of 
Haldane's pseudopotentials^ 



(n,m|yi j2 |n,?7i) 



(11) 



These eigenvalues are summarized in the third column of 
Table |H We can see that the interaction partially lifts the 
degeneracy in the relative angular momentum quantum 
number m within one Landau level. 

To determine the parity of the 4-spinor under particle 
interchange, we perform the following operation^ 



/^AA(z R ,z r y 
I ^ba(zr,Zt 

V&BB(ZR,Z r ); 



^ AA (z R , —Zr 
$Ba(zr, -Z r ) 
^AB(ZR, -Z r ) 

\$bb(zr, -z r )j 



(12) 



where A, B denote the sublattice degree of freedom. We 
thus see that the parity of the two-Dirac-electron eigen- 
states depends not only on the total relative angular mo- 
mentum (m+n) , which fixes the exponent in the variables 
z ri z r in the wavefunction (z r ,z r \n, m) r , but also on the 
second and third components of the 4-spinor. From Ta- 
ble I, we find that, for \^M,m{N = 2) ) with even m, the 
states (I), (III), (VI), (VII) and (VIII) are symmetric un- 
der particle exchange; whereas the states (II), (IV) and 
(V) are antisymmetric under particle exchange. 



B. Wavefunctions resulting by addition of an extra 
particle to the one-particle state 

In addition to calculating the eigenstates and the eigen- 
values E of the two-particle interacting problem, in order 
to compute the DOS in Eq. (1), we also need to construct 
the wavefunction resulting when a particle is added to the 
single-particle state |\&o(-W = 1)): while taking into ac- 
count the overall symmetry. Since the particle-addition 
process is assumed to be instantaneous, and not to per- 
turb the host state, the wavefunction resulting by the 



addition of one particle to the one-particle state can be 
constructed by taking the product of two single-particle 
states and (anti-)symmetrizing it with the symmetriza- 
tion operator Vs or antisymmctrization operator Vas- 
Notice that, here, we need to take into account both 
the symmetrization and the antisymmctrization of the 
orbital wavefunctions. This is because of the spin-valley 
component that can also be antisymmetric or symmetric, 
such that the total wavefunction satisfies fcrmion statis- 
tics. 

Close to v = —2, we take the single-particle state to 
be the n = LL wavefunction (0, |0,mi)), see Eq. (5). 
The addition of an extra electron in the same LL results 
in the two-body state 



I* 



S(AS)/ 



= V. 



S(AS) 





|0,TOll 





|o,m 2 ; 



1 

V2 



(13) 



\|0,mi)i|0,m 2 ) 2 ±|0,m 2 )i|0,mi) 2 / 



where the subscripts 1 and 2 in \n, m)i\n', m') 2 denote 
the subspaces for the respective particles. 

For a generalization to any integer filling n, it is justi- 
fiable to take the (n — 1) LLs to be inert. The two-body 
state I ^5(^45)) can then be constructed in a similar man- 
ner as for the n = case. For example, take the n = 1 
LL case just above the filling v = 2, the instantaneous 
addition of an electron to the n = 1 LL wavefunction 
(|0,mi), |1, ?7ii))/\/2~ results in 



I* 



n— 1 \ 
S(AS)I 



1 

2^S{AS) 



|0,TOi 
l.TOi 



0, m 2 ; 

1, m 2 ; 



(14) 



/|0,mi)i|0,m 2 ) 2 ± |0,m 2 ) 1 |0,m 1 ) 2 N 
1 I |0,mi)i|l,?n 2 ) 2 ± |0,m 2 )ijl,mi) 2 
2y/2 1 II, "H)i|0, 7712)2 ± |l,m 2 )i|0,mi) 2 

\|l,TOl)l|l,TO 2 ) 2 ± |l,rn 2 )i|l,mi) 2 , 

For a translationally invariant system, the DOS calcu- 
lation does not depend on the angular momentum m 2 of 
the added particle, such that we may set m 2 = O^i In 
terms of ladder operators, the two wavefunctions Eq. (13) 
and Eq. (14) are then given by 



l*f(A%) 



\[(b\) m > ± (bir* 



|0,0)i|0,0) s 



(15) 



(b\r- ± (ply 

(b\) mi ± (bh n 



\a\a 



(b\) mi ± (bl) mi 



[{b\) mi ± (b\) mi 



|0,0)i|0,0) 2 , 



(16) 



respectively. By substituting the CM and relative coordi- 
nates operators ai. 2 = (an±a r ) /V2, 6i )2 = (bft±b r )/\/2, 
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and using the binomial expansion, the wavefunctions in 
the new basis become 



\nTAs } ) 



/ y i,mi 



S(AS) 



1=0 



/ y l,mi 



S(AS) 



1=0 



|0,/)ji|0 I mi-0r,(17) 



1 



^(°Jt-°r) 

7f(4 + 4) 

U((4) 2 -(4) 2 )/ 



x |0,0i?j0,mi - l) r , 



(18) 



with the coefficients 

P S(AS) 



^r , -P±(-ir-1 1 /^!(7). («) 

We see that the product state constructed in this way is a 
superposition of the two-body eigenstates of the relevant 
energy level listed in Table H 



C. Total wavefunction 

Having obtained the orbital part of the wavefunc- 
tion for two Dirac particles, we now take into account 
the remaining internal degrees of freedom relevant for 
graphene. The first one is the intrinsic 1/2 spin of the 
electron. The second is the valley index K, K' associated 
to the two inequivalent Dirac cones in the graphene en- 
ergy spectrum. These are separated by a large reciprocal 
wavevector, and even though Coulomb interaction can in 
principle induce inter- valley scattering, the matrix ele- 
ment involving atomic-scale high-momentum exchanges 
are typically small.— Scattering between different valleys 
can thus be neglected in the low-energy regime, such that 
the valley index may be described by a pseudo-spin-1/2, 
with respect to which the Coulomb interaction is approx- 
imately SU(2)-symmetric. The total wavefunction |<&) is 
then a direct product of three parts: 



|$) = |#) <g> | spin, valley). 



from Eq. (13), while we fix r«2 = by invoking transla- 
tional invariance, the groundstate |\&o(Af = 1)) remains 
macroscopically degenerate in the quantum number mi. 
Therefore, the local DOS needs to be averaged over the 
-/Vs-fold degeneracy, e.g., 

, N B -i 

<W=wz E £%-M 



mi=0 a 



|<*a(JV = 2)1^^,1*0^ = 1)>| , (21) 

for the electron-addition part of the DOS. 

In the presence of a high magnetic field, the N = 1 
groundstate is a spin-polarized state with a Zeeman en- 
ergy —A z /2, where A z is the energy splitting between the 
majority and minority spin states. However, the ground- 
state electron can belong to either the K or K' valley. 
Now, when an extra electron is injected, the latter can 
also have a spin pointing either parallel or anti-parallel 
to the ground state electron, and it can reside on either 
the K or K' valley. 

We first consider the groundstate electron residing on 
the K valley and the added electron being spin par- 
allel (S z = 1, where S z is the total spin component 
along the magnetic field direction) but belonging to ci- 
ther of the valleys. It then follows that the total wave- 
function with an added particle can be either |$) = 

<L 1 ,t,K^o,t,K(N = l)) that is 

l*AS°) ® | spin-triplet, valley-triplet), (22) 
or |$) = c i mi yp K ,\^o,t,K(N = 1) ) that is given by either 
|^5 =0 ) <8> | spin-triplet, valley-singlet), 



or 



^AlS ) ® | spin-triplet, valley-triplet). 



(23) 



Here, | ^s(as) ) are given by Eq. (17). Substituting them 
into Eq. (21), and summing over all N — 2 interacting 
eigenstates \^ a (N = 2) ) (partly summarized in Table I), 
we obtain 



(20) 



<y(u-v£=° + A,/2) 



The parity of the orbital part | "J) has been discussed in 
Sect. Ill A and IIIB, whereas the spin and valley parts 
can separately form either a singlet or a triplet state, 
respectively. The total two-body wavefunction must be 
antisymmetric under particle exchange. 



IV. THE DOS OF GRAPHENE FOR n = 

In this section, we compute the electron-addition DOS 
close to v ~ in the n = LL, that is just above the 
filling v = —2 (or the electron-removal DOS close to 
v = 2 by particle-hole symmetry). We first note that 



4- E %-KT° + A z /2) 



N 



(24) 



m£odd 



This is the result which one also obtains in the case of 
2D electrons in GaAs heterostructures, that is the weight 
of the peaks corresponding to odd pscudopotentials is 
3 times larger than that for even ones. However, as 
we shall see, the four-component structure of graphene 
LLs yields eventually a different result than the two- 
component structure in non-relativistic LLs, when S z = 
two-particle states are taken into account. 

In principle, even though the cigenstatc summation is 
performed over all two-body eigenstates, the only one 
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yielding a non-trivial contribution is the eigcnstate (I) of 
Table I. We also note that there is an extra Zeeman en- 
ergy cost of — A z associated with the interacting eigen- 
state because of the spin-triplet component. Further- 
more, to write down the above expression we have em- 
ployed the summation formula 



N B -i 

E 

mi— 



2 mi (mi — m)\m\ 



(25) 



in the thermodynamic limit Nb — > oo, for any integer 
m.— A parallel analysis can be made for the ground state 
electron residing on the K' valley, which leads to the same 
result. 

On the other hand, for the addition of an electron with 
opposite spin, the resulting state with an identical valley 
|$) = cl n i K \^ ,t,K{N = 1) ) gives rise to either 



1*31°) 



| spin-singlet, valley-triplet), 
or 

| spin-triplet , valley-triplet) ; 



(26) 



and the resulting state with an opposite valley |4>) = 
c mi j, K'\^o,t,K(N = 1) ) gives rise to either of the states: 



^AS ) ® | spin-singlet, valley-singlet), 
| spin-triplet, valley-singlet), 
| spin-singlet, valley-triplet), 



l*S =0 > 



I*as°) ® | spin-triplet, valley-triplet). (27) 
The resulting DOS contribution is 
4 



.4 



h=0 i 



N B 



J2 S(w-VZT° -A./2) 



m£even 



n=0 



A./2) 



(28) 



m£odd 



When putting together the contributions from adding 
a parallel-spin electron and an opposite-spin electron to 
the DOS 

4=V) = ^.M+A*$ t=1 {uj), (29) 

the total weight for the odd m peak becomes 5/3 times 
larger than that for the even m (while being possible to 
neglect the Zeeman energy difference). 

As we have already mentioned above, the relative 
weight 5/3 ~ 1.67 between the spectral weight of the 



odd pseudopotentials with respect to the even ones is a 
benchmark of the underlying four-component structure 
of graphene LLs, due to the spin- valley degeneracy. In 
the case of a two-component system, such as in a conven- 
tional 2DEG in GaAs heterostructures, the ratio would 
be 3<2i This result is retrieved in the n = graphene LL 
close to the charge-neutrality point [y ~ 0), where one 
of the spin-valley components (say the spin component 
in the case of a dominant Zeeman effect) is completely 
frozen. The relative spectral weight between the odd 
and even pseudopotentials therefore yields insight into 
the multi-component structure of LLs. 



V. THE DOS OF GRAPHENE FOR GENERAL n 

It is now straightforward to generalize our study to 
other LLs. For a filling factor v ~ 2, we start by consid- 
ering a groundstate which is fully occupied for all n < 1 
LLs and a single spin-polarized electron at n = 1. The 
addition of an electron results in the eigenstate described 
in Eq. (18). Following the same procedure to compute 
the DOS as in the previous section, we obtain 

A%% =1 (oj) = ^- £ <y(w-Oi,m + A,/2) 



m^even 



(30) 



jn^odd 



and 



A%% =q (oj) = — J2 <y(w-Oi,m-A,/2) 



rn^even 



+ ]Vs E 5 ( w - fil -m - Az/2) 



(31) 



m£odd 



where now, only the eigenstate (VII) contributes to the 
DOS, andf2i, m = 2 v / 2fo/! B + (5/8)V':= ,) + (l/4)^ 1 + 
(1/8) V^ =2 . Compared to the usual 2DEG system, we 
see that the position of the peak does not only contain 
information about the n = LL pseudopotcntial V£ =0 
but also on higher LL pseudopotentials V™ =1 ' 2 . This is 
due to the fact that a general two-body eigenstate of the 
interacting problem consists of spinorial components that 
occupy at the same time different LLs n in the relative 
coordinate subspace. 

Let us also write down the solution for the DOS of 
electron addition close to v ~ 6. The two-body state for 
anii = 2 LL electron with an added particle is given by 



7 



I* 



S(AS)/ 



mi 

1 F S(AS) 
o / t l,mi 



1=0 



(at) 2 ) 

K(4) 3 + (4) 3 -(4) 2 4-4(4) 2 ) 

|((4) 4 + (4) 4 -2(a^) 2 (at) 2 ) 

I 



|0,l) H |0,mi-I)T 



(32) 



and taking the overlap with the interacting eigenstatc 
(VIII) from Table U the electron-addition parts of the 
DOS are given by 



41=iM = ]v- E * 

m£euen 

6 



A./2) 



— V S(u - 2 ,m + A,/2) 



(33) 



and 



■ E * 

m£euen 



Az/2) 



£ S(u}-n 2 , m -A z /2). 



(34) 



rn^odd 



where 2 , m = 4fe/Z B + (13/32) + (1/16) V£ =1 + 
(1/4) V™= 2 + (3/16) V^= 3 + (3/32) V™= 4 - Thus, the DOS 
at this filling factor contains rich information on the pseu- 
dopotentials V£ belonging to many LLs. 



VI. DISCUSSIONS AND CONCLUSIONS 

We have calculated the tunneling DOS in graphene in 
high magnetic fields when the filling factor is close to 
v n ± 2(2n +1). In order to describe the electronic in- 
teractions, we have used the method of Haldane's pseu- 
dopotentials, which describes the two-particle interact- 
ing eigenstates of a system in strong magnetic field. The 
method is valid for a system in the very close proximity 
of a completely filled LL, such that, besides an integer 
number of filled and inert LLs, only a single electron or 
hole are present. Although this limit may seem, at first 
sight, extremely theoretical, it describes the experimen- 
tal situation of a very sparsely electron- or hole-filled LL, 
in which the average distance between the particles is 
larger than the cyclotron radius Rq = lBy/2n + S n ,o in 
graphene, that is v -C l/(2n + d n _o). The tunneling from 
the STM tip, which injects a second particle into the sys- 
tem, can thus measure the overlap between the resulting 
state and the two-particle interacting eigenfunctions of 
the system. Also, it allows one to measure the difference 
in energy between the two-particle interacting states of 



the system and the one-particle state, thus yielding in- 
formation about the strength of the interactions. 

Our calculations revealed that the DOS spectrum ex- 
hibit peaks, the energy of which can be related directly 
to the energies of Haldane's pseudopotentials. While the 
n = state is quite similar to the n = LL in non- 
relativistic 2D electron systems with a parabolic band 
dispersion, the higher LL DOS structures are different, 
in that the energies in the spectrum do not result from a 
single pseudopotential value, but involve combinations 
of these values, corresponding to states with different 
angular momenta. This is a direct consequence of two 
graphene-specific properties. The first one is that the 
spinorial eigenstates have (sublattice) components in dif- 
ferent non-relativistic LLs. The second one is that, as 
a consequence of the Lorentz invariance of the underly- 
ing Dirac equation, the center of mass and the relative 
degrees of freedom are intimately coupled. Finally, the 
relative spectral weight between the peaks corresponding 
to odd and even pseudopotentials yields insight into the 
multi-component structure of graphene LLs. 

It would be interesting to generalize our results for 
larger partial fillings, eventually moving towards the 
regime of the fractional quantum Hall effect. Such an 
analysis would allow one to make predictions about the 
experimental spectroscopic signature of different highly- 
delicate quantum Hall states, such as the v — 1/2 state, 
that do not have distinct anomalous transport signatures; 
understanding the nature of such states has been a long- 
standing question in the study of the FQHE. While the 
fractional quantum Hall effect has been measured only 
recently in graphene in transport experiments, spectro- 
scopic measurements may yield additional information 
about the relevant electronic interactions in graphene 
LLs and thus about the nature of strongly-correlated 
phases in partially filled levels. 
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TABLE I: Two-body eigenvalues and eigenstates 
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